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Abstract 



The oxygen deficient CuO x planes of the high-T c superconductors 
REBa2Cu30g+x (RE denotes a rare-earth atom or yttrium) are studied. 
Starting from the Emery model for the Cu-0 chain fragments, an extended 
ASYNNNI model is proposed to describe the statistical mechanics of the 
CuO x planes. The model is analysed by using Monte Carlo techniques, 
computing especially the charge transfer to the Cu02 planes regulating the 
high-T c superconductivity, and the concentration of differently coordinated 
Cu ions. Results are compared to experimental data for various RE cases. 
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Typeset using REVI^X 



I. INTRODUCTION 



The high-T c superconductors REBa 2 Cu 3 6+x (REBCO), where RE denotes yttrium 
or one of the rare-earth atoms Yb, Er, Y, Sm, Nd, or La, have been studied extensively, 
both theoretically and experimentally. Among others, the impact of the oxygen deficient 
CuO x (sometimes called Cu(l)) planes on superconductivity has been investigated. In 
particular, the charge transfer from the Cu(l) planes to the superconducting Cu0 2 (or 
Cu(2)) planes has been shown to play a crucial role. 

The Cu(l) planes consist of linear Cu-0 — • -O-Cu chain fragments (CFs). Each 
Cu atom is strongly coupled to the two apical oxygen atoms in adjacent Ba-0 planes. 
Depending on the number, I, of neighboring oxygen atoms in a CF, Cu atoms may be 
(/ + 2)-fold coordinated. Obviously, a threefold coordinated Cu atom, Cu^ 3 \ is located at 
the end of a chain fragment; fourfold coordination (Cu^) is realized for a copper atom 
in the CF, while an isolated Cu atom is twofold coordinated, Cu 1 - 2 ). The concentration of 
differently coordinated Cu atoms has been determined experimentally by means of NMR- 
NQR Hm and x-ray absorption spectroscopy (XAS) Using the soft XAS method, it 

is possible to estimate the concentration of oxygen holes in sufficiently long Cu-0 chains. 
For example, for untwinned YBa 2 Cu 3 7 , it has been found to be about 30% ||. 

T c , the superconducting phase transition temperature, as a function of oxygen content 
x is known to display the prominent 60 K and 90 K plateaus for RE = Yb, Y and Er 
0, i.e., the compounds of comparatively small RE ionic radii. The analog of the 60 K 
plateau is less pronounced for RE = Gd J7|, and the plateaus are practically absent for 
the REBCO compounds with RE = Nd |],[7j and La ||. The distribution of chain lengths 
in the CF ensemble allows one to estimate the concentration of holes leaving the chains, 
rih, which is usually assumed to be proportional to T c , unless one is in the underdoping 
(n/ l <0.05) and overdoping (n^>0.25) regimes. Underdoping corresponds to localization 
of transferred holes in Cu0 2 planes, while overdoping holds near the maximum of T c (nh) 
(see @ and |0|). 

T c depends not only on x but also on the temperature, T q , at which the Cu(l) planes 
have been effectively equilibrated experimentally. T q is typically the annealing temper- 
ature for rapid quenches, or the room temperature, Tr, if the cooling was done slowly; 
below T R , the mobility of the oxygen atoms in the CuO x planes is greatly reduced. Vary- 
ing T q up to 350°C in experiments on YBCO, Veal and co-workers [O] showed that the 60 
K plateau is sensitive to changes of T q . In fact, the plateau becomes less pronounced with 
increasing T q . A subsequent experiment [I2~| confirmed these results, and demonstrated 
that the various annealing temperatures T q are accompanied by rearrangements in the 
ensemble of chain fragments. Note, that the temperature considered in our theoretical 
analysis of the Cu-0 chains corresponds to T q in experiments. 

The Cu-0 chains in the Cu(l) planes may form differently ordered phases, the or- 
thorhombic, Ortho-I, Ortho-II, and probably Ortho-Ill (see, for instance, |T3|-p^1, as well 
as the tetragonal phases, which have been identified experimentally. The perfect Ortho- 
I and Ortho-II structures together with a typical arrangement of oxygen atoms in the 
tetragonal structure are shown in Figure |l[ 

The transformation from the semiconducting to the metal state in YBCO is usually 
ascribed to the tetragonal-orthorhombic transition. The underlying microscopic mecha- 
nism seems to be the derealization of holes in the Cu0 2 planes resulting in a change of 
the oxygen hole chemical potential, which, in turn, influences the structural properties 
Pi. 
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Some years ago, a lattice-gas model, the asymmetric next-nearest-neighbor Ising 
(ASYNNNI) model, has been proposed for describing the oxygen ordering in the Cu(l) 
planes of YBCO [19|]. It is based on three effective coupling constants as shown in Fig- 
ure 0. Two of these couplings, the nearest neighbor oxygen repulsion V\ (> 0) and the 
next-nearest neighbor attraction V2 (< 0) mediated by a copper atom, favor Ortho-I-like 
configurations. The third coupling, V3, the direct interaction between next-nearest oxy- 
gen atoms without intermediate Cu atom, is supposed to be repulsive, thereby favoring 
formation of Ortho-II-like configurations. The phase diagram of the model is sketched in 
Figure |^. As noted before, when comparing to experiments, the temperature should be 
interpreted as T q . 

The phase diagram of the ASYNNNI model has been analyzed using various numerical 
approaches, such as cluster variation methods [^-22], Monte Carlo simulations for static 
|T7| . |23| -p5f as well as dynamic properties ||26|| , and transfer-matrix finite-size techniques 



|27j. Among the approximate analytical approaches, one may mention a low-temperature 
expansion of the free energy J28| and high-temperature expansions |29| , |3"0|j . The model 
allows to reproduce remarkably well the main features of the experimental phase diagram 
for YBCO. As a note of caution, attention may be drawn, however, to experimental find- 
ings providing evidence for a possible discrepancy in the location of the phase boundary 



of the Ortho-II phase |T2| , |3T| -p^] . According to these experiments, the disordering of the 
Ortho-II phase, oxygen content x close to 0.5, occurs at about 450 K. 

In the ASYNNNI model the internal degrees of freedom of the CFs are ignored, such 
as the charge and spin distributions. To take them into account, one may start from an 
appropriate microscopic model of the Cu-0 chains and try to map it onto a corresponding 
lattice-gas model, possibly of ASYNNNI-type. (Microscopic models have been also em- 
ployed by Aligia @,[3^] to explain and predict some superstructures in the Cu(l) planes). 
Such an approach has been attempted before JS^rrSj, motivated by the suggestion that, 
in the case of YBCO, the ASYNNNI model may be insufficient to describe phenomena 
in the part of the phase diagram where short Cu-0 chains dominate. For instance, the 
chemical potential isotherms 0] are reproduced only poorly, whose singularities seem to 
be important near the tetragonal-to-Ortho-I phase boundary where the average lenght 
of the chains does not exceed a few lattice constants. It had been suggested p8|-4IO|] to 
describe the tetragonal phase in YBCO at x < 0.4 in terms of a phase separation. In 
particular, neutron and x-ray data had been interpreted, in the regime of phase separa- 
tion, as signalling short chains with a herringbone structure. However, subsequently, the 
superstructure reflections were attributed to a parasitic phase (Nevertheless, there 
may be a chance to find herringbone structures or some other regular configurations of the 
very short chain fragments in LaBCO, but not in YBCO; for possible patterns, see Figure 
|j. Indeed, it has been shown [41] that those chain fragments are electric quadrupoles, and 
if they are dominant, then the herringbone structure would be energetically preferred). 

The motivation of the present study is to reconsider the ASYNNNI model in view of 
the new experimental facts. We shall relate that lattice-gas model to a strongly correlated 
electron model, the Emery model, for the Cu-0 chains. This will allow us to reinterpret 
the crucial coupling V2, which will turn out to depend, for instance, on temperature. 
In addition, a new parameter will be introduced, reflecting the relevance of the shortest 
chain fragments, i.e. Cu-O-Cu. Using Monte Carlo simulations, properties of the result- 
ing extended ASYNNNI model will be computed, based on the statistical ensemble of 
chain fragments. The main considerations in relating the Emery model to the ASYNNNI 
model will be presented in the next section, followed by a section in which, based on 
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the lattice-gas model, some interesting properties like the temperature dependence of the 
charge transfer from the Cu-0 chains in the Cu(l) planes to the Cu(2) planes for various 
REBCO compounds as well as the distribution of differently coordinated Cu atoms will 
be discussed. A short summary will conclude the article. 



II. MICROSCOPIC DESCRIPTION 

Charge transfer mechanism 



Following prior considerations ||42|| , we shall describe briefly a possible charge transfer 
mechanism for the Cu-0 chains in the Cu(l) planes. 

Assuming a strong repulsion between neighboring oxygen sites, see Figure [|, linear Cu- 
O chain fragments are formed with essentially no crossings. In a CF, oxygen and copper 
atoms form a strongly correlated system due to a strong hybridization of the oxygen p x (p y ) 
orbitals with the d z 2_ x i{d z 2_ y 2) orbitals of the copper atoms. If we ignored the oxidation 
process, the n oxygen atoms in a CF of composition Cu n+ iO n would be neutral, while the 
n+1 Cu atoms would be in the monovalent state, as in the parent compound YBa2Cu306- 
In reality, a charge transfer is expected to occur, which may be described formally as a 
two-step process. The first step corresponds to a charge redistribution in the CF, leading 
to only one divalent 2 ~ ion in each CF; schematically 
Cu+-0-...-Cu + -0-Cu+-...-0-Cu + — ► 
Cu 2+ - 0-- . . . - Cu 2 +- O 2 -- Cu 2+ - . . . O Cu 2 +. 

The second step involves a charge transfer from the monovalent oxygen ions in a chain 
fragment to the Cu(2) planes or to other fragments. The resulting CF has the form 
Cu 2 + 1 2 n "0~_ m , where m — 1 holes have been transferred, assuming that each copper ion 
remains in the divalent state. It has been shown experimentally || that for long chains in 
YBCO m/ra~0.7. For short chains, the optimal CF configurations have been estimated 



H to be Cu^+O^O- for n = 4, Cu^+Oj^Cr forn = 3 and Cu| + 2 -Q- (or Cu^O 2 ," 



for n = 2. Accordingly, no hole (one hole) leaves the n = 2 CF, one hole and two holes 
leave the n = 3 and n = 4 CF, respectively. Generalizing this finding into a simple rule of 
thumb, one may estimate m to be, for an arbitrarily long chain, m = nint (n(m/n) opt ), 
where nint denotes the nearest integer. (m/n) opt is that optimal ratio, which is realized 
in long CFs. 

The Emery model 

As the starting point for a microscopic description of the Cu-0 chains we use the 
one-dimensional Emery model [43], applied to systems of finite length. The Hamiltonian 
is 

H = -t \Pla d P° + (fyaPro) + Yl ^V U r + U p n r fn rl ) + 

<r,p> a r 

+ (^n p + U d n p] n pl ) + U cnd J2 n p (!) 
p p 

where the sum runs over the O sites (r) and the Cu sites (p) of a finite chain, < r, p> 
denotes neighboring lattice sites, and p refers to the Cu® sites at the chain ends. d a and 
p a [d\ and p a ) are annihilation (creation) operators of a hole with spin a on Cu and O 
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sites, respectively. n a is the hole (spin a) occupation number and n = Y, a n a- This model 
ignores the contribution of the apical oxygen ions to the energy, which, indeed, play a 



minor role, because they are essentially divalent (cf discussion in pH}|). 

Because of the quantum origin of alternating Cu-0 chains, the hole occupation number 
on Cu sites is not rigorously equal to 1, as tacitly assumed when the charge transfer 
mechanism was discussed. This occupation number fluctuates, remaining slightly below 1 
(~ 0.85, as numerically estimated in Now, a good quantum number is the number 

of holes (both, oxygen and copper), say v{n) in a CF of length n. Nevertheless, m remains 
a quantum number, characterizing the hole transfer from a chain. The relation between 
m(n) and v{n) is given by 

u(n) = In- (m(n) - 1). (2) 



There are strongly correlated models, e.g., the Kondo- reduced Emery model (see p2| , p6| ). 
which allow to treat Cu ions of CFs as divalent, and oxygen ions as monovalent (n — m) 
or divalent (m). 

Note that we include in model ([[]) the shift of the on-site hole energy level, U en d, which 
is due to non-equivalency of Cu*- 3 - 1 and Cu^ 4 - 1 ions. There may be alternate ways to take 



into account the U en d-teim; for instance Aligia uses a model [0,0] which includes the 
interaction of charges of neighboring atoms, p(Cu) and r(O), 

U pd £ (l + 72 p )(-2 + n r ) 

<r,p> 

Our description could be interpreted as a simplified version of the C/ p d-term, U en d ~ 
£7 p d(2 — where ?v denotes the thermal expectation value for the occupancy of an 
oxygen site. 

In the case of RE=Nd or La, the ionic radii of the rare-earth ions are larger than in 
the other REBCO compounds. Consequently, the interatomic distance is also larger, but 
the hole hopping amplitude is smaller - in the chains, and in the planes. We shall discuss 
below why this circumstance may favor short chain fragments, depending on the oxygen 
hole chemical potential, /i, which varies in the REBCO series. 



Free energy 

We consider the Cu^ 1 0^"0~_ m chain fragment, described by the Emery model ([[]). 
We may determine, say, numerically, its free energy (here and in the following, the free en- 
ergies are taken per oxygen atom), f(n; m). Indeed, performing statistical averaging with 
respect to the internal, fermionic, degrees of freedom (m and possible spin arrangements), 
and disregarding the configurational entropy, the free energy for long chains (n ^> 1) has 
been recently shown to have the form [|3(] 



4>(n) « 0i Inn. (3) 

n 2n 

where 0i and 02 are coefficients which are given by the energy parameters of the Emery 
model; 02 is, approximately, linear in T. 



In specifying the parameters of the Emery model, we largely follow Refs. |45],[T8|. 
is taken as the reference level. Measuring the other energy parameters in units of t, we 
set e p = 2, Up = 4 and Ud = 8, similar to the values in the Cu02 planes of YBCO. The 
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REBCO compounds are expected to differ mainly in the value of the hopping amplitude 
t (being approximately t = 1.35 eV « 15000 K for YBCO). As noted before (@^[[8|, the 
remaining two parameters, the chemical potential \i and U cn d, play an essential role in 
favoring long or short chains, fi equilibrates the hole exchange between the two oxygen 
hole reservoirs, CuC>2 planes and CFs. U en d stems from the non-equivalency of the Cu^ 3 -* 
and Cu^ 4 ) sites (ed + U en( ± denotes the hole level on a Cu^-site). Although f/ end is positive, 
as seen from its relation to £/ p d, its exact value is not known, neither from experiments 
nor from band-structure calculations. 

In YBCO the oxygen atoms of the oxygen deficient planes tend to arrange themselves 
into long chains at low enough quenching temperatures; probably even at quite low oxygen 
content. In contrast, in LaBCO long chains are predominant only at x ~ 1 0. In the 
REBCO compounds, the various tendencies in the chain lenghts are expected to be caused 
by the various ionic radii of their rare-earth constituents (see Table |). In the model 
description ([I]) of those compounds, the crucial parameters are the hopping amplitude t 
which depends on the RE ionic radius, and the chemical potential which regulates the 
degree of filling of the ID O-hole band. 

To see whether short or long Cu-0 chains are favored, we compute the free energy 
per oxygen atom as a function of the chain length. If it increases with length, then short 
chains are more likely, otherwise longer chains are preferred. 

For a given set of parameters one may calculate the free energy in the grand canonical 
ensemble, thus incorporating the chemical potential /i, which governs the concentration 
of holes in a chain. Following this scheme one needs, in principle, the whole set of energy 
levels of the Emery model. However, since the temperatures we are dealing with are 
much lower than the typical electronic energies, which are of order several thousand K, 
it is sufficient to include the lowest energy levels in the free energy calculations. These 
levels may be determined by using the Lanczos algorithm j|7|, as we did in evaluating the 
free energy of chains with up to n = 6 oxygen atoms (plus 7 copper atoms). The detailed 
procedure for calculating the free energies is described in Appendix A. 



III. THE EXTENDED ASYNNNI MODEL 

Estimate of effective parameters 

The free energy of the chains described by the Emery model, expression (^), corre- 
sponds in the standard ASYNNNI model to the energy of a Cu-0 chain with n oxygen 
atoms, neglecting interchain interactions, i.e. 

4>{n) =V 2 - — . (4) 
n 

A clear difference of 4>(n) and 4>{n) is that <fi(n) decreases monotonically with n at V 2 < 0, 
whereas, if the analogous condition 02 < is satisfied, <p(n) exhibits a minimum at finite, 
but exponentially large value n opt oc exp20 2 /7V It is estimated as exceeding a hundred 
at temperatures below 300° C. Hence, the logarithmic correction in the r.h.s. of Eq.(|3]) 
may start playing an important role only at x ~ 1. 

02 in Eq.(|3|) depends on temperature, while V 2 does not (in both cases, we did not 
consider the configurational entropy). Accordingly, accepting the Emery model as the 
more fundamental description for the chain fragments, the ASYNNNI model may be 
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modified by introducing an effective, temperature-dependent coupling V 2} corresponding 
to 02 in the Emery model. The difference in the coefficients <\>\ and 02 may not break 
the analogy with the ASYNNNI model. In relevant experiments, the CF ensemble is 
measured at a fixed oxygen content, thus 0i contributes equally to any realization of CFs 
and may be rescaled, say, to the value of 02- 

Note that Eq.(^) holds only for sufficiently long chains. Indeed, for short chain frag- 
ments, 0(n) may show oscillations fl8f . To include these subtle features in a statistical 
description (for example, in a Monte Carlo simulation), one may deal with the free en- 
ergies of the short chains separately. The largest deviation from the expressions (3) and 
(4) occurs at n = 1. Thence, to incorporate those oscillations in the simplest fashion, we 
consider an extended ASYNNNI model, attributing to the chains with n = 1, Cu-O-Cu, 
an additional energy V 2 . For larger values of n, we keep the parameters of the standard 
ASYNNNI model. Accordingly, the free energy may be cast in the form 

$( n ) = V 2 -—, n>l; i(l) = V 2 . (5) 
n 

To substantiate Eqs.(Q) and (|5|), we shall now present results of calculations of the free 
energy for the Emery model of Cu-0 chain fragments. 

Let us first roughly estimate a realistic value of U en d- In Figure the free energy 
of the Emery model is shown, using units appropriate for YBCO (we have chosen the 
chemical potential \x at each U cn d to assure that the concentration of oxygen holes in 
the longest chains we studied is ~ 0.3). Specifically, the energy parameters are taken in 
units of t = 1.35 eV ~ 15000 K, so that the temperature range of experimental interest 
(T q between 300 K and 1200 K) is between 0.02 and 0.08. Comparing the free energies 
for long chains to expression Eq. (|J) (disregarding, at this stage, V 2 ) one finds that 
U en d — 1 corresponds to a moderately attractive V 2 , of several hundreds of Kelvin in 
agreement with the experimentally observed disordering of the Ortho-II phase at about 
450 K. Otherwise, for example, for f/ end = the free energy grows too rapidly with the 
chain length, corresponding to an effective repulsion V 2 of the order of a few thousands 
of Kelvin (see also |3f|). Increasing U en d leads first to a weaker repulsion, and finally 
an attractive V 2 . For instance, U en d = 1.5 would imply a value for V 2 of the order of 
several thousand K, being a strong effective attraction. Thereby, by tuning U en d, one may 
stabilize either short chain fragments (some realizations at x < 1/2 are shown in Figure 
H) or very long Cu-0 chains. 

While U en d is expected to be essentially a constant for all REBCO compounds, the 
chemical potential [i may vary, e.g., with the oxygen content x or the RE element. In 
Figure |5] the free energy per oxygen is shown for several values of the chemical potential at 
the particular choice U en d = 1. With \i increasing, the effective intra-chain interaction V 2 
between the oxygen atoms in the Cu-0 chains changes from attraction to repulsion. This is 
happening in a fairly broad range, [1=1.5 - 2.0, where V 2 goes from a moderately attractive 
value (~ —0.15) to a moderately repulsive one (~ 0.1). Changing the RE element in the 
REBCO compounds implies changing [i and hence the intra-chain interaction. As will be 
discussed later, V 2 varies slowly, remaining ~ —0.1 in the range of relevant values of [i. 
In fact, V 2 introduces an important competing effect. 

In Figure [7] we depict the free energy per oxygen atom versus n for U en d = 1 at different 
temperatures, choosing the chemical potential so that the concentration of oxygen holes is, 
for long chains, approximately 0.3, as it is the case for YBCO. /i ~ 1.6 is well compatible 
with that concentration. One may easily see that the lower the (quenching) temperature 
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T the larger the effective intra-chain oxygen-oxygen attraction will be. In addition, at low 
T, the free energy shows an oscillatory behavior for small values of n (see, for example, 
Figure [7| at T = 0.02). Such oscillations are maximal at T = (see also Figure 1 in ||18|| ). 
It is obvious from Figure [7] that the fitting of the free energy to that of the standard 
ASYNNNI model, Eq.@, is unsatisfactory even in the high temperature region. The 
n — 1 case has to be fitted separately, defining V 2 in the extended ASYNNNI model, 
Eq.(||). It should be emphasized, that the effective parameters of the hereby defined 
extended ASYNNNI model, V 2 and V2, depend not only on the temperature, but also on 
U en d- In addition, they depend on a specific filling of the ID oxygen hole band, i.e. \i. 

The other two couplings of the full ASYNNNI model, V\ and V3, are of predominantly 
electrostatic origin, and may be supposed to be constant. Of course, their values do not 
follow from the Emery model. (At this point, attention may be drawn to another attempt 
to determine parameters of the (standard) ASYNNNI model from first-principles, based, 
in that case, on linear muffin-tin-orbital calcualtions for YBCO | 48fl . In that analysis, 



only ground state properties have been considered, leading to constant couplings V\, V 2 , 
and V3.) 

In Figure [7|, also least-square fits of the free energy for the Emery model to the ansatz, 
Eq. (5), defining the extended ASYNNNI model are shown, at several temperatures. 
From such fits one may determine V2 and V2 as function of temperature. Examples are 
shown in Figures |8| and |9] for a few values of U en d around U cn d = 1 , presuming a fixed hole 
concentration of 0.3 in the longest chains. It is evident from Figure || that V2 displays a 
pronounced linear dependence on T (cf, Eq.([3]), where 02 is also linear in T), and V 2 is 
practically linear in U en d at a given temperature. 

To fix t/end) one may compare results on the extended ASYNNNI model to experiments. 
A possible procedure will be demonstrated now for YBCO. There we estimate U en d by 
locating the disordering transition of the Ortho-II phase of YF^CusOg+x with x = 0.5, 
known to occur experimentally at about 450 K |T2]j3l| -|33l ■ That transition temperature 
is found to depend rather sensitively on the choice of U en d, and hence V% and V2 in the 
framework of the extended ASYNNNI model. The lattice gas model is analysed by using 
standard Monte Carlo (MC) techniques. 

To simulate the extended ASYNNNI model, one first has to specify the couplings V\ 
and V3. Different values have been proposed and studied for the conventional ASYNNNI 
model |Hj23,22]. Here we adopt values which have been argued to reproduce the ex- 



perimental structural phase diagram of YBCO fairly well [||],[^|, namely, in terms of 
V2, 

Vi = -2.86 ■ V 2 (6) 

V 3 = -0.46 • V 2 (7) 

As discussed above, V 2 is expected to depend on temperature much more strongly than Vi 
and V3. Thence, in fixing V\ and V3, we choose the value of Vj> at the Ortho-II disordering 
temperature, Tn ~ 0.03 (note that we checked, by simulations, that choosing the room 
temperature Tr, instead of Tn, leads to rather minor changes of V\ and V3, of a few 
percent, and affects the estimate for U en d only mildly). 

Simulations at T n for different values of U en d (and tuning /1 to reproduce the hole 
concentration for the longest chains), which determine V^Tn; U en d) and ^(T n ; U en d) have 



been performed. The order parameter ran [23,p5| of the Ortho-II phase is found to drop 



abruptly at U cn d = 1.02, see Figure |10|. From Figures |8| and |9|, one obtains V^Tn) ~ —0.12 
and V 2 (Tu) rs —0.09. Summarizing, we have finally estimated, for YBCO, U en d = 1.02 
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and \i = 1.59 from the experimental values of Tfi and the hole concentration for long 
chains. 

In substituting yttrium by another RE element, the chemical potential will be changed. 
This, in turn, will lead to different values of the effective parameters in the ASYNNNI 



model. Results of pertinent calculations are depicted in Figure |TT], taking U en d = 1.02. 
By changing fi in between 1.55 and 1.70, one induces only mild changes in V2, in between 
-0.15 and -0.1. More importantly, we observe in that range a crossover in the tendency to 
form short or long Cu-0 chains. Obviously, that tendency follows from the fact whether 
V2 is smaller than V2, or vice versa. At \x = 1.59 (characteristic for YBCO), the long 
chains are favored; at \x ~ 1.65, a compensation holds. At larger values of the chemical 
potential, short chains will be preferred. 

Accordingly, for each REBCO compound, one may predict the preference for short or 
long chains from the value of /i. Although the chemical potential is not precisely known 
for the various RE cases, one may safely assume that it increases with the ionic radius. 
As shown in table |, these radii vary appreciably. For instance, they are larger for Nd and 
La than for Y. 

A microscopic description for estimating \x should invoke the structure of the oxygen 
hole bands of the chain fragments (ID) as well as of the Cu(2) planes (2D). It has been 
assumed in [O] that the ID oxygen hole band is lower than the 2D band, with an overlap, 



as shown schematically in Figure [12|, where /1 is above the bottom of the 2D hole band. 



Despite a potentially significant charge transfer from chains (i.e. m(n) — 1 from each 
CF of length n), such a band structure is no contradiction. In fact, at small oxygen 
content x, shorter CFs are preferred, the charge transfer gets effectively surpressed. The 
holes transferred from chains become partly localized on the impurity levels and partly 
2D carriers, with a rather low, semiconducting, concentration. In the metallic state the 
chemical potential intersects the 2D band; a sufficiently large hole transfer from chains to 
planes causes an electrostatic shift of the levels, and the ID band raises with respect to 
the 2D band. Based on such a description, the transformation from the semiconducting 
to the metallic state, by changing the oxygen content x, can be qualitatively described 
0®- 



In principle, one could estimate the variation of \x with x from a self-consistent scheme 
which would be very tedious, based upon experimental data which have not been es- 
tablished reliably, so far. In the following analysis, we shall, however, assume that the 
chemical potential is independent of x. This will be sufficient to reproduce the main 
trends resulting from the change in fi in the REBCO series. To describe the correspond- 
ing thermal properties of the chain fragments, we shall investigate the extended ASYNNNI 
model, with the parameters V2 and V% changing with the chemical potential. For V\ and 
V3, we shall always take the values for YBCO, see above. Because they are associated 
with Coulomb forces on small distances, they are expected to vary only slightly, by a few 
percents, when substituting Y by, say, Nd or La. As we checked in simulations, such a 
variation will not affect the conclusions. The results of the Monte Carlo study will be 
compared to experimental findings for LaBCO and NdBCO |{],@|. 

Monte Carlo simulations, Comparison with experiments, Discussions 

The simulations of the extended ASYNNNI model describing REBCO compounds 
have been performed for systems with L x L Cu sites, imposing full periodic boundary 
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conditions. Usually, L was taken to be 40; to estimate the role of finite size effects, a 
few runs were done with larger lattices. We employed Kawasaki dynamics, keeping the 
oxygen content x constant during a simulation, and choosing the initial and final site of 
an oxygen atom randomly. The new site, if empty, is accepted with a probability given 
by the Boltzmann factor of the change in energy of that possible move. An updating 
scheme which restricts the hopping of the oxygen atoms to neighboring lattice sites would 
be more realistic in mimicing dynamical processes, but slower. Because we are here only 
interested in equilibrium properties, we employed the faster algorithm. 

By varying the oxygen content x and the temperature T, the structural phase diagram 



may be determined. Figure [13] depicts the phase diagram for the case of YBCO, using the 
values of the effective parameters discussed above. The phase boundaries have been esti- 
mated from the behavior of the order parameters of the various phases and the positions 



of the maxima in the specific heat, as before [p3[| . The topology of the phase diagram is 



similar to that obtained for the standard ASYNNNI model in prior computations, e.g. 
23| , |2"4|| , including the Ortho-I, Ortho-II, and tetragonal phases. Perhaps most notably, 



the maximal phase transition temperature of the Ortho-II phase is, however, somewhat 
reduced, in accordance with the experimental findings. 

By analysing Monte Carlo equilibrium configurations one may easily identify the Cu-0 
chain fragments and hence determine, for instance, the thermally averaged concentration 
of differently coordinated Cu ions. In Figure [14] the concentration of Cu*- 2 -* ions, C2, is 



shown as a function of the oxygen content for various temperatures, in the YBCO case. 
While at high temperatures C2 decreases monotonically with increasing oxygen content, 
one finds a non-monotonic or plateau-like behavior when lowering the temperature, so 
that one may encounter the Ortho-II phase. Obviously, in the perfect Ortho-II structure, 
at x = 0.5 and sufficiently low temperatures, half of the Cu atoms belong to the perfect 
Cu-0 chains extending throughout the lattice and half of them are in oxygen-free rows, 
see Figure ]l](c). Accordingly, the concentration of Cu^ 2 -* ions in such an ideal configuration 
is 0.5. Experimental results of XAS measurements at low and high temperatures M are 
included in Figure |14[ Agreement between experimental data and theoretical curves is 
quite satisfactory, except for the point at x ~ 0.39, low T q , which may, indeed, have 
been determined inaccurately in the experiment. The behavior of c<i versus x at T = 450 
K also agrees rather well with the experimental results using NMR-NQR fl]||. The 
concentrations of Cu^ 2 \ Cu^ and Cu^ 4 ^ (c 2 , c 3 , and C4) obtained from the MC simulations 
at T = 450 K in the case of YBCO are shown in Figure [H| together with experimental 
data (their accuracy has not been mentioned) 0. A few experimental points deviate 
clearly from the theoretical curves. Note, that the experimental data of plotted in 
Figure [15] obey well the sum rule 

c 2 + c 3 + c 4 = 1 (8) 

This identity assumes that all Cu ions in the Cu(l) planes can be identified as Cu*- 2 -*, Cu^ 3 -* 
or Cu 1 - 4 -*. With that assumption, a second sum rule is evident 

c 3 + 2c 4 = 2x. (9) 

These two sum rules allow to express two of the concentrations, say, c 3 and C4, by the 
remaining one, C2, and the oxygen content x; for instance, 

c 3 = 2 — 2x — 2c2, C4 = 2x — 1 + C2. 
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The concentrations are allowed to vary in the range defined by the following inequalities, 

0<x<l/2 — ► 1 -2x < c 2 < 1 -x, 0<c 3 <2x, < c 4 < x; . . 

1/2 < a; < 1 — ► < c 2 < 1 - x, < c 3 < 2 - 2x, 2x - 1 < c 4 < x. ^ ' 

A few points of the experimental set disagree significantly with (|9|), several points are 
beyond the limits defined by (|T0D. Note that there is an uncertainty, reaching up to 0.1 
according to ||, in determining the oxygen content in the range 0.3 <x <0.4 (probably, 
in the region of phase separation). 

As mentioned above two concentrations can be determined by the remaining third 
concentration. Taking the latter one as c 2 for x < 0.5, c 3 for 0.5 < x < 0.7, c 4 for 
0.7 <x < 0.8, and again c 2 for 0.8 <x (c 4 would be beyond the limit of the inequality), 
we can fit nicely the revised set of experimental points to theoretical curves in Figure |15|. 

Similar measurements on the concentrations of differently coordinated Cu ions have 
been reported for NdBa 2 Cu 3 6+x (NdBCO) and LaBa 2 Cu 3 6+x (LaBCO) @. As argued 
above, we propose that the effect of the increase of the ionic radius in Nd and La com- 
pounds can be incorporated into the microscopic model by an increase in the chemical 
potential, leading to modified effective parameters in the extended ASYNNNI model. Al- 
though /i is expected to depend on the oxygen content x for each REBCO compound, 
we consider a constant fi for a given RE, for reasons stated before. This simplification 
seems to be justified in order to detect the main trends. Indeed, when increasing /i from 
1.59, for YBCO, to 1.62 and 1.66, for NdBCO and LaBCO, respectively, the simulated 



data on c 2 , c 3 , and c 4 agree rather well with the experimental results, see figures 16 and 



17| . For NdBCO we can use c 3 for determining the other concentrations in the range 
0.5 < x <0.7, as in the YBCO-case. For LaBCO, c 3 and c 4 have not been accessible for 
direct measurements in the range 0.5 < x; in this case c 2 plays the determining role. 

There may be also five- and sixfold coordinated Cu ions in Cu(l) planes, due to a 
crossing of perpendicular chain fragments. For example, Cu*- 5 -* would appear, if an oxygen 
atom formed a 90° configuration with two other oxygen atoms (all three are supposed 
to be adjacent to the same Cu atom). The energy loss due to a Cu^ (Cu^) ion would 
be 2Vi (4Vi). From the simulations, the concentration of such ions is seen to be very 
small for all oxygen contents x. A kink in a CF would be energetically somewhat more 
favorable, albeit being still rather costly, V\. At the kink, a fourfold coordinated Cu ion 
occurs, although it differs from the conventional Cu^ 4 -* ion. Such a kink configuration 
cannot easily be described by the Emery model. 

Another important quantity related to the chain fragments of the oxygen deficient 
Cu(l) planes, is the charge (hole) transfer, rih, to the Cu0 2 planes. It depends on the 
chemical potential as well as on other external parameters such as x and T q . ^(n), the 
expectation number of holes (thermal averaging) leaving a CF of length n, should be 



m(n) — 1 (see Figure [18], showing a rather non-trivial dependence on n); it is associated 



with the expectation number of oxygen and copper holes i/(n) through (fj). For long CFs 
in YBCO, one has m(n)/n « 0.7; but for short chain fragments, it may differ significantly. 
The average concentration of oxygen holes transferred from the chains, the charge transfer 



n h , may be approximated 



1 JL m(n) - 1 



n h (fi;x;T q ) = — J2 5 <A^ C f(^)>, (11) 



n=l 



where < Acf(^) > is the average number of CFs of length n. The factor of 1/2 appears, 
because holes from one Cu(l) plane are supplied to two Cu(2) planes. 
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Note that one can calculate the number of holes v{n) in the framework of the Emery 
model, at least for fairly short chains, see Appendix A, which may then be extrapolated 
to longer chains. 

Figure [19] shows the charge transfer versus oxygen content for a set of equidistant 
temperatures, obtained from simulating the extended ASYNNNI model for the YBCO 
case. One observes that below x w 0.25 the charge transfer to Cu02 planes is extremely 
small. In this regime, in the tetragonal phase, short chain fragments dominate, contribut- 
ing much less to the charge transfer than the long chains. The preference for short chains 
could be reduced by going to very low (quenching) temperatures which, however, are 
usually not accessible experimentally (recall that T q is bounded by Tr from below). The 
tetragonal phase is still semiconducting, and the holes are mainly localized in Cu02 planes 
in the vicinity of the CF from which they originated. The formation of longer CFs leads 
to a significant increase of the hole transfer, giving rise to the semiconducting-to-metal 
transformation. 

The charge transfer displays a plateau-like behavior near x = 0.5, reflecting the well- 
known 60 K plateau in the transition temperature to the superconducting phase in YBCO 
at that oxygen content. This aspect has been discussed in quite a few publications be- 
fore. In many experiments the close relationship between the superconducting transition 
temperature T c and rih has been found |50|J5T[] . For instance, T c was argued to be linear 
|52] , [53|| or quadratic in rih |54]j55|]. In both cases, the plateau in T c would correspond to 
one in n^. A second, so-called 90 K plateau in T c versus x occurs at x > 0.8, where rih 
increases monotonically, while T c remains almost constant. Overdoping is believed to be 
responsible for that plateau PJ54}|: T c practically saturates when n h becomes larger than 
the optimal hole concentration (« 0.2 — 0.25). 

In Figure ^ the charge tranfer rih versus oxygen content x is depicted at T = 450 K for 
three values of the chemical potential, mimicing, as before, YBCO, NdBCO, and LaBCO, 
i.e. \x = 1.59, 1.62 and 1.66. One readily sees, in accordance with the experiments, that 
the low-temperature plateau in T c versus x is expected to be much less pronounced in 
NdBCO [[F]|| and absent in LaBCO ||. In that figure, we also included Monte Carlo data 
for ii = 1.56, mimicing a REBCO compound where the rare earth ions would be even 
smaller than yttrium, like YbBCO (see table |). Indeed, it was shown experimentally 
that the 60 K plateau in YbBCO is even more manifest than in YBCO, in accordance 
with the MC data. 

MC simulations allow to compute easily the average length of Cu-0 chain fragments, 
I. Results are displayed in Figure Ell. t is given by the ratio (total length of CFs) / (total 



number of CFs). The numerator is identical with the number of oxygen ions, while 
the denominator is, neglecting kinks and crossings of CFs, one half of the Cu® atoms. 
Thence, one would obtain I = 2x/c^. That formula has been used in JJ, to express the 
experimental findings for C3 for YBCO, NdBCO, and LaBCO in terms of £. However, care 
is needed in doing that. In Figure JH], we compare Monte Carlo results for I invoking the 
approximate formula to those using the correct expression for the average chain lenght. 
Indeed, the overall agreement is fair, but there are obvious differences near x ~ 1. 



IV. SUMMARY 

The present paper has been motivated, first, by recent experimental data on Cu-0 
chain fragments 0, and, second, by the wish of elaborating and concretizing the idea to 
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relate a microscopic description of those fragments to a statistical analysis, i.e. to relate 
the quantum mechanics of individual chains to the statistical mechanics of their ensembles 
fl8|j . Our main considerations and findings may be summarized as follows. 

• We consider the Emery model, which belongs to the family of strongly correlated 
fermionic (hole) models, incorporating into it the effect of chain ends ([/end-term). 
We assume that the lowest, Cu-like, band is completely filled, whereas the interme- 
diate, O-like, band is partly occupied. The occupancy of the latter band is regulated 
by the chemical potential fi, which determines the equilibrium of the oxygen hole 
reservoirs in the Cu02 planes and CFs. 

• In accordance with experimental results, the chemical potential level is assumed to 
be located in the ID oxygen hole band of the Cu-0 chain fragments. The oxygen 
hole concentration in long CFs is determined by the chemical potential fi. 

• The oxygen hole number in a CF follows from the number of holes leaving the chain. 
There is essentially no charge transfer from short CFs of lengths n = 1 and 2. This 
explains the predominance of fairly short CFs in the semiconducting state of the 
REBCO compounds. 

• The "semiconductor- metal" transformation is governed by the dependence of the 
chemical potential on the oxygen content x. Analysis of finite chains by means 
of the Emery model shows that CFs of minimal length occur less often when /i 
decreases. This is expected to be accompanied by an increasing number of holes 
transferred from chains to planes. Coulomb forces lead to the lowering of fi, which 
favors longer CFs, giving rise to an additional charge transfer, etc. 

• In our analysis, the REBCO compounds have been characterised by different, but 
fixed values of fi. However, in reality, with x increasing, /i is believed to decrease 



proportionally to n h , i.e. the hole doping of the Cu0 2 planes. Using Figure [H] , 
these changes of /i could be incorporated into the calculations. 

An important role in the Emery model is played by U eri <i] the energy loss stemming 
from the ends of a CF determines the effective attraction or repulsion between the 
oxygen atoms in CFs. In fact, ignoring this term would result in a very strong 0-0 
repulsion in the chains (of the order of several hundreds meV in the case of YBCO, 
where the concentration of oxygen holes in long CFs is about 0.3). Experimental 
data allow one to estimate U en d to be ~ 1 in units of the hopping term t of the 
Emery model. 

The free energy of the Emery model has been calculated for individual CFs, char- 
acterized by their length function of the chemical potential /i, Z7 en d, and 
the temperature T. For specific REBCO compounds, /i and U en d can be fixed us- 
ing experimental data, e.g., for the oxygen hole concentration and the disordering 
temperature of the Ortho-II phase. By comparison of the free energy, the values 
of the parameters of an extended ASYNNNI model can then be estimated. In the 
framework of that model, the structural phase diagrams of the REBCO compounds 
can be conveniently studied, using, for instance, Monte Carlo simulations. 

The larger the ionic radius of the REBCO compound the higher the chemical po- 
tential, and the shorter the Cu-0 chain fragments. 
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• Analysing the extended ASYNNNI model for various RE constituents, we also com- 
puted macroscopic properties such as the concentration of differently coordinated 
Cu ions and the charge transfer from the chain fragments. 

In conclusion, the theoretical approach presented in this paper (start from microscopic 
calculations, identify effective parameters of a statistical model as well as, based on exper- 
imental results, their values, and finally perform simulations), allows to reproduce main 
trends in the statistical properties of Cu-0 chain fragments throughout the various RE- 
BCO compounds. Refinements and extensions, based on additional experimental data, 
are possible. 
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APPENDIX A: THE FREE ENERGY CALCULATION FOR CHAIN 

FRAGMENTS 

The Hamiltonian of the Emery model for a CF, characterized by n and m, deals with 
v = 2n + 1 — m fermions which are spin-1/2 holes. Positive integers m (> 1) describe the 
charge transfer from a CF. The eigenstates of this Hamiltonian can be easily divided into 
blocks which correspond to different numbers u [y = 0, 1, . . . , 2(2n+l)) and, at given u, to 
different total spin projections (S z = —v/2, —v/2 + 1, . . . , v/2). Any block can be further 
subdivided into two blocks which consist of either symmetric or antisymmetric states with 
respect to the central Cu atom of a CF. The size of actual blocks grows fastly with n: 
CFs with n = 6 oxygen atoms, plus 7 copper atoms, are still accessible for calculations of 
a substantial number of the low lying energy levels. 

Our conjecture is that the CF properties, including the concentration of holes, depend 
on the chemical potential // which, in turn, is a function of the hole doping. Thus, we 
must calculate the grand canonical partition function: 

2(2n+l) v/2 

2n(T,/i)= £ J2 5>-[*C.**>~H/r (Al) 

v=0 S z =-u/2 i 

where % runs over all the states within the block at given v and S z . The temperatures of 
interest do not exceed 1200 K, which is low on the energy scale of the model (|l|). Indeed, 
e.g., including one thousand low lying energy levels at n — 6, we can properly evaluate 
the partition function (doubling the number of low lying states, one does not change the 
free energy significantly). 

According to |J, one has rn/n ~ 0.7, provided the copper hole occupation number 
is equal to one. We incorporate this fact, demanding the average total hole number 
V = 7 + 6 • 0.3 ~ 8.8 for the longest chains (n = 6 oxygen and seven copper atoms), 



analysed using the Lanczos algorithm [47]. Thus, v is defined as 
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1 2(2n+l) i//2 

17 =^ E E 5>e-^^-^/ fcT (A2) 

i.=0 S z =-v/2 i 

and /i(T) is then numerically determined from the constraint V = 8.8 for n = 6. 

In the framework of the Lanczos algorithm M71, applied to a block of the states {u,S z }, 



we can find, say, the M lowest eigenvalues. AM iterations should guarantee convergence 
of M eigenfunctions. There may be a problem with the Lanczos algorithm: it fails to 
distinguish two or more degenerate eigenvalues. However, this is not a serious problem 
in our case, since the accidental degeneracy occurs only rarely in the Emery model. This 
has been checked by calculating exactly the eigenvalues for CFs up to n = 3. 
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TABLES 



Atom 


Ion 


Z 


Radius 


Ytterbium 


Yb 3 + 


70 


0.86 


Yttrium 


y3+ 


39 


0.89 


Neodymium 


Nd 3+ 


60 


0.995 


Lanthanum 


La 3 + 


57 


1.02 



TABLE I. Ionic radius in A and atomic number Z of some trivalent rare earth ions RE 3+ . 
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FIGURES 



FIG. 1. Shown are two possible domains of the Ortho-I structure with orientations of chains 
along the horizontal (a) and vertical (b) axes; one of four domains of the Ortho-II structure (c); 
and a typical pattern of the tetragonal configuration with orientations of finite chain fragments 
along both x and y axes (d). Squares and filled (open) circles denote Cu atoms and O atoms 
(vacancies) , respectively. 



FIG. 2. Sketch of the ASYNNNI model, with interactions coupling oxygen atoms. Squares 
and filled (open) circles symbolize copper and oxygen atoms (vacancies). 



FIG. 3. Schematic structural phase diagram of YBCO. Below the dashed curve, i.e. below 
Tfj, equilibrium can hardly be reached experimentally. 



FIG. 4. Examples of simple perfect structures consisting of the shortest CFs at x = 0.5, (a) 
and (b). The herringbone structure (a) would be favored by pure electrostatic (quadrupolar) 
interactions [41|. Squares and filled (open) circles denote copper and oxygen atoms (vacancies). 



FIG. 5. Free energy per oxygen atom, /, of Cu-0 chains versus chain lenght (or number of 
oygen atoms), n, for various values of t/ en d at T = 0.06 ~ 900 K. For each curve, the chemical 
potential has been chosen so that the concentration of oxygen holes is 0.3 for the longest chains, 
independent of f7 e nd- The free energy is measured in units of the hopping parameter t 15000 
K. 



FIG. 6. Free energy per oxygen atom, /, of Cu-0 chains versus length, n, for equidistant 
values of the chemical potential, ranging from 1.5 to 2.0, at U en d = 1.0 and T = 900 K. 



FIG. 7. Free energy per oxygen atom, /, of Cu-0 chains vs length, n, for a few temperatures, 
with U en d = 1.0 and the hole concentration being 0.3 for the longest chains. Open symbols denote 
fits to expression (|5|). 



FIG. 8. Vi, vs temperature, T, for a set of equidistant values of U en d- The hole oxygen hole 
concentration is taken to be 0.3 for the longest, n = 6, chains. 



FIG. 9. V2 vs temperature, T, for a set of equidistant values of U en ^. The oxygen hole 
concentration is 0.3 for the longest chains. 



FIG. 10. Ortho-II order parameter, mn, vs U en d, at T\\ = 0.03 corresponding to 450 K. The 
full symbols denote Monte Carlo data for the extended ASYNNNI model, applied to YBCO; 
the interpolating line is included as a guide to the eye. 



FIG. 11. V2 (full symbols) and V2 (open symbols) vs chemical potential [i for various tem- 
peratures, at [/end = 1-02. 
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FIG. 12. Illustration of shifting of bands by going from the semiconducting (tetragonal) to 
the metallic (orthorhombic) state in REBCO compounds. The ID hole band filling is larger in 
the semiconducting state. 

FIG. 13. Phase diagram (temperature T(K) in Kelvin vs. oxygen content x) obtained from 
MC simulations of the extended ASYNNNI model for YBCO. Typically, systems with 40 x 40 
copper ions were considered. 

FIG. 14. Concentration of twofold coordinated Cu ions, C2, versus oxygen content, x, at 
various temperatures. Simulations of the extended ASYNNNI model for YBCO were done at 
T = 300 K and 600 K. MC data are shown by full dots; interpolating lines are included, for 
clarity, as guides to the eye. Experimental results are denoted by open symbols. 

FIG. 15. Concentration of two-, three- and fourfold coordinated Cu ions, q (I =2,3, and 4) 
versus oxygen content, x, at T = 450 K, simulating the extended ASYNNNI model for YBCO, 
(j, = 1.59 (full symbols and, as guide to the eye, lines). For comparison, experimental data |l]] 
are shown (open symbols). 

FIG. 16. Notation as in Figure 15; however, in this case for NdBCO, taking the chemical 
potential to be 1.62. 

FIG. 17. Notation as in Figure 15; however, in this case for LaBCO, taking the chemical 
potential ijl to be 1.66. 



FIG. 18. Number of holes per oxygen atom, i/^(n)/n, leaving a CF of length n (n = 2, . . . , 6), 
as function of temperature T, at fj, = 1.59 (y^ = for n = 1). 

FIG. 19. Charge transfer, n^, as function of oxygen content x, at various temperatures from 
300 K to 600 K, simulating the extended ASYNNNI model for YBCO. 

FIG. 20. Charge transfer, rih, as function of oxygen content, x, simulating the extended 
ASYNNNI model mimicing RE= La, Nd, Y and Yb (from bottom to top). 

FIG. 21. Average chain length, £, as function of oxygen content, x, simulating the extended 
ASYNNNI model mimicing RE= La, Nd, and Y (from bottom to top). Full and open symbols 
denote direct and indirect (through C3) determination of £, respectively; see text. 
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